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ABSTRACT 

Analysis of 20 months of observations at the Parkes radio telescope shows secular changes in the 
pulsed emission from J0737— 3039B, the 2.77 s pulsar of the double-pulsar system. Pulse profiles are 
becoming single-peaked in both bright phases of the orbital modulation although there is no clear 
variation in overall pulse width. The shape of the orbital modulation is also varying systematically, 
with both bright phases shrinking in longitude by ~ 7° yr . However, the combined span of the two 
bright phases is relatively constant and together they are shifting to higher longitudes at a rate of 
~ 3° yr^ 1 . We discuss the possible contributions of geodetic precession and periastron advance to the 
observed variations. 

Subject headings: binaries: general — pulsars: general — pulsars: individual (PSR J0737— 3039B) 
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1. INTRODUCTION 

One of the many uniq ue features of the double-pulsar 
syste m J0737-3039A/B (jBurgav et al.ll2uM ILTne et all 
|2004|) is the dramatic orbital modulation of the pulsed 
emission from PSR J0737-3039B (hereafter "B"), the 
2.77 s pulsar of the system. Both the shape and the in- 
tensity of the B pulse profile vary with orbital longitude, 
with two bright phases centered around longitudes (with 
respect to the ascending node) of ~ 210° (bright phase 1, 
hereafter bpl) and ~ 280° (hereafter bp2), respectively. 
The pulse is weakly visible at other orbital longitudes, 
except perhaps between 30° and 60°. 

This unprecedented behavior has been interpreted as 
being due to the interaction between radiation from the 
23 ms pulsar J0737-3039A (hereafter "A") and B's mag- 
netosphere (I.Tenet fc Ransomll2004l iZhang fc Loebll2004l 
lLvutikovll2004U2005lL Direct evidence of the mutual in- 
teraction between the two pulsars is manifested by the 
dri fting behavior at A's pe riod in the single pulses of 
B l)McLaughlin et al1l20fl4al) and by the modulation at 
B's period of the emissi on from A during its eclipse 
l|McLaughlin et alJl2004b^ . 

The geometry of the system is expected to vary 
with time as a result of the strong relativistic ef- 
fects occurring in the system: the geodetic precession 
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ijDamour fc Ruffinl 119741) and periastron advance of B 
have periods of only 71 and 21 yr, respectively, both 
muc h less than for an y other known pulsar binary sys- 
tem ijLvne et al.ll2004|) . These effects are expected to re- 
sult in_se£ula£_changes in the observed pulse shape (see 
e.g. lKramerlll9 98). in the interaction between the two 
pulsars, and hence in the orbital modulation. Somewhat 
surprisingly, Parkes observations of A show no significant 
evidence of profile-shap e variation over a 3 yr interval 
({Manchester et aLll2005|L Near alignment of A's rota- 
tion axis with the orbit normal is a possible explanation 
for the observed lack of variations. 

In this Letter we report the detection of secular varia- 
tions in both B's pulse shape and the orbital modulation 
of B's emission using observations with the Parkes ra- 
dio telescope spanning 20 months. In 321 we present the 
observations and results and in []3]we discuss their inter- 
pretation. 

2. OBSERVATIONS AND RESULTS 

The double-pulsar system J0737— 3039A/B was ob- 
served at 20 cm (1390 MHz) using the Parkes 64 m radio 
telescope from 2003 May 1 to 2005 January 1. All obser- 
vations were made using the central beam of the multi- 
beam receiver, except from November 2003 to September 
2004 when the "H-OH" receiver was used. A filter-bank 
system providing 512 total-power channels across 256 
MHz was us ed. Further details of the observing systems 
are given in [Manchest er et alJ l)2005j) . For the present 
analysis we used a total of 94 hr of data spread over 105 
observations, more or less uniformly distributed across 
the 20 month data span. 

Figure ^ shows the intensity of the two bright phases 
of emission from B as a function of the orbital longi- 
tude with respect to the ascending node for six epochs 
across the 20 month data span. Each panel was ob- 
tained by summing the observations taken over an ^3 
month interval. Th e PSRCHIVE data analysis system 
( Hota n et al.l [2.004) was used to sum and image the 
da ta. These p l ots c onfirm t he overall picture presente d 
in iLvne et"aT1 l|2004fl and iRamachandran et a 
showing the complex pulse shape and flux density vari- 
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Fig. 1. — Intensity of the PSR J0737— 3039B pulse emission at 1390 MHz as a function of orbital longitude and pulse phase at six epochs. 
Only the longitude range 180°-310° covering the two bright phases bpl (lower) and bp2 (upper) and a pulse-phase window of 0.1 pulse 
periods centered on the pulse are shown. The dashed lines represent the longitude of periastron at each epoch. 

ations as a function of orbital longitude. Furthermore, 
Figure Q clearly shows that both the extent in orbital 
longitude of the bright phases and the shape of the pulse 
at a given orbital longitude are evolving with time. We 
also note that the mean pulse intensity in bp2 is decreas- 
ing relative to that in bpl. In the following sections we 
quantify these trends. 

2.1. Pulse Profile Variations 

As FigureQJshows, the mean pulse profile for B is quite 
different in bpl and bp2 and even varies somewhat with 
orbital longitude within each bright phase. Despite the 
small variations within each bright phase, we have cho- 
sen to illustrate the trends in profile shape by integrat- 
ing across bpl and bp2 separately for each of the epochs 
corresponding to those in Figure Q The resulting pulse 
profiles are shown in Figure |2 This figure shows sig- 
nificant long-term variations in the mean pulse profile 
for both bright phases. The variations are in fact very 
similar in both bpl and bp2, with the leading compo- 
nent getting weaker and the trailing component becom- 
ing broader with time. There is little or no significant 
evolution in the 10% width of either profile (upper lim- 
its of 7 and 8 ms yr _1 for bpl and bp2, respectively). 
For bpl, the width of the trailing component appears 
to increase with time, but the change is nonlinear, with 
the first two profiles narrow and the last four broader. 
The feature on the leading edge of the trailing compo- 
nent in the 2004 January profile suggests that this width 
increase is due to the relatively rapid growth of an ad- 
ditional component between the two other components. 
For bp2 the 50% widths are difficult to interpret as the 
leading component has a relative amplitude close to 50% 
and the signal-to-noise ratio is smaller. 

2.2. Orbital Phase Variations 

To quantify the orbital evolution of the bright phases, 
we took the data illustrated in Figure Q and summed the 
pulse intensity at each epoch over the 3% of pulse phase 



Fig. 2. — Mean pulse profiles, normalized in peak, of PSR 
J0737-3039B at 1390 MHz for the two bright phases, bpl (left, 
orbital longitude 190°-235°) and bp2 (right, orbital longitude 260°- 
300°) at six epochs from 2003 June (top) to 2004 November (bot- 
tom). All profiles have been scaled to the same peak height. Num- 
bers near each profile are the pulse widths at 10% (left) and 50% 
(right) of the pulse peak. For bpl, these widths have an uncertainty 
of about 1 ms, and for bp2, about 3 ms. 

where the pulse is present. A baseline level calculated 
in the orbital longitude range 30° - 60° was then sub- 
tracted, and the area under the resulting curves equal- 
ized to remove the effects of refractive interstellar scin- 
tillation. The resulting normalized pulse intensity curves 
are shown in Figure [3] for each of the six epochs. 

This figure shows clear systematic changes in the or- 
bital modulation with time. The most obvious change 
is that the saddle region between the two bright phases 
is becoming less filled: that is, the end of bpl and the 
start of bp2 are becoming more separated. Less obvi- 
ous are that the modulation as a whole is moving toward 
later longitudes and that relative to bpl, bp2 is becoming 
fainter. These changes are shown more quantitatively in 
Figures 0] and [3] 
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Fig. 3. — Variation of pulse intensity as a function of orbital 
longitude for six epochs. The small bar represents the typical rms 
noise in the baseline. 
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Fig. 4. — Variations with time of the start and end of the two 
bright phases, measured at the 10% points, and the centroids of 
the two phases separately and together. The dotted line shows the 
precession of periastron during the 20-month interval. 
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Fig. 5. — Variation with time of flux(bpl/bp2), the relative 
integrated pulse intensity of the two bright phases. 

Linear least-squares fits to the points in Figure _] show 
that the starting point of bpl and the end of bp2 are 
slowly shifting toward higher longitudes at a rate of 
3?4 ± 1?7 and 2? 5 ± 1?7 yr" 1 , respectively (all the er- 
rors are quoted at 2 a). Likewise, the centroid of bpl 
and bp2 taken together is moving at a rate of 3?0 ± 1?2 
yr _1 , confirming that the overall modulation pattern is 
shifting toward later longitudes. However, the interior 
points have a different behavior, with the end of bpl 
moving backward at — 3? 8 ± 1?7 yr _1 and the start of 
bp2 advancing at 9?6 ± 1?7 yr _1 . This shows that bpl 
and bp2 are shrinking at a rate of 7? 2 ± 2?4 and 7? 1 ± 2? 4 
yr _1 , respectively, and the gap between the two bright 



phases is widening at 13? 4 ± 2? 4 yr _1 . 

Figure _] shows the relative energy in bpl and bp2, 
obtained by integrating under the curves in Figure _] 
from 190° to 235° for bpl and from 260° to 300° for 
bp2, and taking the ratio. Clearly, bp2 is getting fainter 
relative to bpl; a linear fit to the ratio has a slope of 
l.OiO.lyr- 1 . 

3. DISCUSSION 

This Letter reports the first observation of secular 
changes in the pulsed emission from the J0737— 3039A/B 
binary system. The B pulse profile, initially consisting 
of two obvious components in both bright phases, has 
become single peaked over the 20-month interval of our 
observations. This constrasts with th e apparent stabil- 
ity of A's profile over a similar interval l)Manchester et all 
2005). Moreover, the orbital extent of both bright phases 
is shrinking at a rate of ~ 7° yr _1 and bp2 is getting 
progressively fainter with respect to bpl. Both the outer 
edges and the centroid of the two bright phases taken 
together are advancing in orbital longitude at a rate of 
~ 3° yr _1 , whereas there is no significant variation in 
the separation of the outer edges of the modulation pat- 
tern. All these numbers have been derived by analyzing 
observations made at 1.4 GHz; similar trends appear in a 
preliminary inspection of unpublished data taken at the 
Green Bank and the Giant Metrewave Radio Telescopes 
at lower frequencies. 

If these changes were to continue at the current rate, 
the second bright phase would disappear by the end 
of 2007 when its centroid will be at orbital longitude 
~ 305°, whereas bpl will fade away in mid-2009, with 
the centroid at orbital longitude ~ 210°. We note that 
the shape of the integrated pulse profile of B at orbital 
longitudes away from the two bright phases is single- 
peaked with 10% and 50% widths (~ 65 and ~ 30 ms, 
respectively) similar to the corresponding values for the 
pulse profile in bpl. This may indicate that B's emission 
is tending to become more uniform along the orbit as the 
physical processes producing the bright regions cease to 
operate. 

Provided B's spin and total angular momentum vec- 
tors are not aligned, gen eral relativity predicts (e.g., 
IBarker fc O'Connell 1197^1 that B's rotational axis will 
precess about the orbit normal at a rate £l g = 5?1 yr _1 . 
This precession will vary the angle of impact of A's 
wind on B's open magnetic field lines. Several authors 
(ILvne et alJ200l Uenet fc RansompOOl ILvutikovll2004t 
iZhang fc Loebll2004HLvutikovll2005|) have suggested that 
the pulsed emission from B is modulated by the imping- 
ing energetic flux from A, so some kind of evolution in 
B's orbital light curve may be expected. We note that 
Q g is of the same order of magnitude as the rate of the 
observed variations in the location and extent of bpl and 
bp2. 

Geodetic precession will also change our line of sight 
across B's emission cone, affecting the shape and width 
of the observed pulse profile as reported for other pulsar s 
fWeisherg et a.lJlTflSfl IKrameHlTlwst IStairs et al.l f20041 . 
In fact, the profile shapes in both bpl and bp2 are clearly 
evolving. However, there is no significant variation in 
the angular extent of the pulse emission, with limits of 
about 1° yr _1 on the change of the 10% widths. This 
may indicate that B's spin axis is almost parallel to 
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orbit normal, possibly supporting the hypothesis that 
wind-torques from A, which dominate the energetics of 
the system, have caused the spin axis of B to almost 
align with the direction of the orbital angul ar momen- 
tum l|Arons et alJl200flt iDemorest et al.ll2004D . 

Periastron precession has moved the orbit with respect 
to the ascending node by 28° over our data span. This 
large angular shift is not reflected in any of the measured 
variations of the orbital light curves and so the relative 
position of the two stars with respect to periastron can- 
not be a significant factor in the changing orbital mod- 
ulation of the emission from B. Likewise, it cannot be 
responsible for the shrinking width of the bright phases. 
In fact, despite the longitudinal extents of bpl and bp2 
undergoing a similar decreasing trend, the linear distance 
between the two pulsars has changed with opposite signs 
in bpl and bp2 along the data span: an ~ 3% increment 
for bpl and an ~ 1% decrement for bp2. 

Various models have been proposed to explain the 
dramatic orbital modu lation of the B pulse emission. 
Uenet fe Ransom! l)2004l) suggested that B brightens only 
when the emission beam from A illuminates B's magneto- 
sphere. In its present form, the model i s inconsistent with 
the o bserved stability of A's emission ijManchester et all 
2005), and it does not include B's geodetic precession or 
the effect of the varying separation of the two pulsars at 
a given orbital longitud e. 

iZhang fc Loebl l)2004|) propose that the modulation is 
a consequence of particles from A's wind streaming into 
some of the open magnetic field lines of B. The model re- 
quires that A's wind is anisotropic and also that B's spin 
axis is not parallel to the orbital angular momentum vec- 
tor, thus allowing geodetic precession to have an effect: 
the size and the location of the portion of B's open mag- 
netosphere directly seen by A's wind at any given orbital 
longitude would be modified. 

Alternatively, iLvutikov l &2W l[j first suggested that B 
brightens when the line of sight runs parallel to the mag- 
netic field lines in the portion of the B's magnetosphcrc 
where the field lines are swept back owing to the im- 
pinging A's wind: in fact, strongly curved magnetic field 
lines would favor the growth of a Cerenkov-drift insta- 
bility, producing radio emission beamed along (and po- 
larized perpendicular to the plane of) the magnetic field 
lines. In this model geodetic precession affects the loca- 
tion along the orbit (and the spin phase of B) at which 
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the condition of alignment between the bent-back field 
lines and the lin e of sight is met. 

More recently, ILvutikovl l)2005|) has proposed that the 
modulations of the B light curve are due to orbital phase- 
dependent distortions of B's magnetosphere by A's wind. 
In this model, B is intrinsically bright at all orbital 
phases, but strong emission is detected only at those 
orbital longitudes where the magnetic field lines along 
which the radio emission is generated are bent toward 
the observer. 

Assuming the geom etry inferred in 

iLvutikov fc Thompson! IpM a ~ 75° and S ~ 60°, 
where a is the angle between spin and magnetic axes 
of B, and S is the angle between spin axis and orbit 
normal) , this model roughly accounts for both the extent 
and the location of the two bright phases. It can also 
accou nt for the variations of the pulse shape alon g the 
orbit ijLvne et al.ll2004t iRamachandran et aLll2005|) and 
predicts that the separation between bpl and bp2 will 
vary because of geodetic precession. The magnitude and 
direction of the changes we observe will place further 
constraints on this model parameters. 

With continued monitoring of the secular evolution 
of the emission from the J0737— 3039A/B double-pulsar 
system, including the measurement of the changes in po- 
larization properties, we should be able to remove the 
existing degeneracy both in the orientation of the spin 
and magnetic axes of the two pulsars and among vari- 
ous precessional effects; this will also be crucial for mea- 
suring the moment of inertia of pulsar A once we mea- 
sure higher-order terms to periastron advance at the sec- 
ond Post-Newtonian level ijKramer eTaOHoog). Com- 
bining these data with an improved understanding of 
the complicated magnetospheric interactions and emis- 
sion physics in the system, we may hope to eventually 
obtain an explanation for this most intriguing pulsar be- 
havior. 
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